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We performed a search for short gravitational wave bursts using about 3 years of data of the 
resonant bar detectors Nautilus and Explorer. Two types of analysis were performed: a search for 
coincidences with a low background of accidentals (0.1 over the entire period), and the calculation 
of upper limits on the rate of gravitational wave bursts. Here we give a detailed account of the 
methodology and we report the results: a null search for coincident events and an upper limit 
that improves over all previous limits from resonant antennas, and is competitive, in the range 
h rss ~ 10~ 19 , with limits from interferometric detectors. Some new methodological features are 
introduced that have proven successful in the upper limits evaluation. 

PACS numbers: 04.80.Nn, 95.30.Sf, 95.85.Sz 



I. INTRODUCTION 



In the quest for gravitational waves (GWs), a primary 
role among the possible sources has always been played 
by those astrophysical events that are expected to pro- 
duce GW bursts, such as the gravitational collapse of 
stars or the final few orbits and the subsequent coales- 
cence of a close binary system of neutron stars or black 
holes. The search for such transient GW requires the 
use of a network of detectors. In fact, the analysis of 
simultaneous data from more detectors at different sites 
allows an efficient rejection of the spurious outliers, cither 
caused by transient local disturbances or by the intrin- 
sic noise of the detectors. Resonant GW detectors have 
operated for decades in several laboratories around the 
world, reliably staying on the air for long periods with 
high duty cycle [l-Q, mainly looking for burst events. 
The coming of age of laser interferometer detectors [H, 
with much better sensitivity and bandwidth, has lead to 
a gradual phasing out of many resonant detectors. 

The ROG Collaboration has built and operated two 
cryogenic, resonant-mass detectors, EXPLORER 0-0| 
at CERN and NAUTILUS @, at the INFN Frascati 
National Labs (Italy). Both detectors have been on the 
air since the early '90s, performing various joint coinci- 
dence searches fUH12| . In the period May 5th 2005 to 
April 15th 2007 they took part in the IGEC2 network 
Q that collected and exchanged data, together with 
the Auriga detector at the INFN Legnaro National Labs 
(Italy) fla ] and with the Allegro detector at LSU (USA) 
jl4| . After that period, Allegro was shut down and data 
have been collected by the three surviving antennas, but 
never analyzed before. 



All these detectors use the same principles of opera- 
tion. The GW excites the odd longitudinal modes of the 
cylindrical bar, which is cooled to cryogenic temperatures 
to reduce the thermal noise and is isolated from seismic 
and acoustic disturbances. Both Explorer and Nautilus 
consist of a large aluminum alloy cylinder (3 m long, 
0.6 m diameter) suspended in vacuum by a cable around 
its central section and cooled to about 2 K by means 
of a supcrfluid helium bath. To record the vibrations 
of the bar first longitudinal mode, an auxiliary mechan- 
ical resonator tuned to the same frequency is bolted on 
one bar end face. This resonator is part of a capacitive 
electro-mechanical transducer that produces an electri- 
cal a.c. current that is proportional to the displacement 
between the secondary resonator and the bar end face. 
Such current is then amplified by means of a dcSQUID 
superconductive device. Nautilus is also equipped with 
a dilution refrigerator that enables operations at 0.1 K, 
further reducing the thermal noise. In the period con- 
sidered, however, the refrigerator was not operational, in 
order to maximize the detector duty cycle. Both detec- 
tors are equipped with cosmic ra y te lescopes, to veto ex- 
citations due to large showers |15Hl7j . The two telescopes 
rely on different technologies (scintillators for Explorer, 
streamer tubes for Nautilus) but both provide a monitor 
of comparable effectiveness and a continuous check of the 
antenna sensitivity (lil Il9j . 

At present, while the large interferometers VIRGO and 
LIGO are undergoing massive overhauls to upgrade their 
sensitivity, there still are two resonant detectors, Nautilus 
and Auriga, that continue to operate in " astro- watch" 
mode, i.e. as sentinels recording data that could be ana- 
lyzed in conjunction with a significant astrophysical trig- 
ger, such as the explosion of a nearby supernova, or any 
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FIG. 1: Spectral sensitivity curves of EXPLORER and NAU- 
TILUS. The two bandwidths overlap for a large fraction of the 
total sensitive region. 



lectcd by an automatic event finder procedure and con- 
stitute the database for our analysis. 

The paper is organized as follows: in sect|TT] we de- 
scribe the data collected in these 3 years, the filtering 
procedure, the criteria chosen to segment the total ob- 
servation period in 5 subperiods and the vetos applied 
to both data and events. In sect JIIII the procedures of 
software injection and time-delayed coincidences are de- 
tailed. The detectors are characterized in terms of effi- 
ciency and accidentals in the various time segments and 
the ROCs are generated, both for the individual anten- 
nas and for the combined observatory; on the basis of 
these parameters, in sect II VI the thresholds for the "true" 
(zero-delay, on time) coincidence search are chosen and 
the search is performed. Finally, in sect|Vl we describe 
the procedure used to compute the upper limit for the 
flux of GW radiation. This procedure, quite different 
from those used in the past, is optimized in each of the 
subperiods and for each of the amplitudes of GW signals 
considered. 

Some final considerations conclude the paper. 



astronomical event thought to be a possible source of 
GW. 

We report here a study on three years of data from Ex- 
plorer and Nautilus, starting from the end of the IGEC2 
network, April 16,2007 and stretching till June 10, 2010, 
when Explorer ceased operations. The spectral sensitiv- 
ity of the two detectors is shown in fig HI 

The purpose of this paper is to describe a search for 
short burst coincident events in the 3 years of data. The 
main interest of this analysis lies in several novelties 
that were implemented in the data analysis procedure 
and that are here detailed; namely, the construction of 
.Receiver Operating Characteristics for the 2-detector ob- 
servatory, in addition to the ones for each antenna, the 
optimization of the threshold pairs in order to maximize 
the detection efficiency, given an a-priori choice of the 
background of accidental coincidences and an optimized 
procedure for the calculations of upper limits on the rate 
of GW burst. The search was carried out keeping a low 
level of accidental coincidences, that we set at 0.1 over 
the entire period. Along with this search, we also per- 
formed a calculation of the upper limit (UL) on the rate 
of delta-like GW pulses impinging on the Earth. The 
method here described presents relevant improvements 
with respect to previously published searches performed 
with resonant detectors: we have used software injec- 
tions of known signals to measure the efficiency of detec- 
tion for each antenna and for the combined observatory. 
Based on the efficiencies so evaluated, and on the mea- 
sured rate of accidentals, the analysis parameters were 
separately optimized for the coincidence search and for 
the UL evaluation. 

Throughout this paper, we shall call "events" or "out- 
liers" those data points in the filtered data stream that 
are larger than a given threshold: these points are se- 



ll. THE DATA 

Data are collected by the two detectors with almost 
identical hardware and software. The output of the 
SQUID amplifier is conditioned by band pass filtering 
and by an anti-aliasing low-pass filter, then sampled at 5 
kHz and stored on disk. Sampling is triggered by a GPS 
disciplined rubidium oscillator, also providing the time 
stamp for the acquired data. 

The data are processed off-line, applying adaptive, fre- 
quency domain filters. We first "whiten" the data, i.e. 
remove the effect of the detector transfer function, pro- 
ducing the so called reconstructed h. A filter matched 
to delta (or very short) excitations is then applied to 
this stream. We shall call hs(t) the output of this pro- 
cess, to remind that the time series so produced does not 
represent a generic wave amplitude h, but is the best 
effort to detect an impulsive excitation. The noise char- 
acteristics estimate is updated averaging hg(t) over 10 
minutes periods. The filters used in the above procedure 
are computed using a model for the detector, fitted with 
the measured values of frequencies and decay times of the 
system resonances, and the experimental noise spectrum. 
Both the model and the signal response were validated by 
hardware injections of known signals: the filtered output 
matched the expected value to better than 10%. 

As usually done for resonant detectors @, the filtered 
output h$(t) is normalized assuming that the excitation 
lasts ~ 1 ms and has a bandwidth of ~ 1 kHz centered 
and flat in the region of the detectors sensitivity. If these 
assumptions are fulfilled, the peak value of hg(t) gives 
the h rss of the input signal. 

The filter is designed and optimized for delta-like sig- 
nals, but it works equally well [2(| for a wider class of 
short bursts, like e.g. damped sinusoids with decay time 
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r < 5 ms. Typical GW signals of this kind and their 
possible astrophysical sources have been discussed and 
exposed e.g. in [2(| and references therein. 

Although the detectors produce quite stationary data, 
their characteristics did change a few times over such 
a long observation period: in some instances, these dif- 
ferences were due to actual changes in hardware (e.g. 
substitution of a preamplifier), other times to some non 
identified factors. We found it useful, to the purpose of 
the study detailed below, to segment the analysis in dif- 
ferent periods where both detectors had noise behavior 
(average noise energy) [24| consistently stable, within the 
statistical fluctuations. This allows us to better optimize 
the search in each period. Consequently, we identified 5 
time stretches (see table I), that roughly coincide with 
solar years, and ran separate optimized analysis on each 
subperiod. Stretch #2, the end of 2007, covers a short 
period, when Nautilus operations were badly disturbed: 
we shall show that the adopted procedures automatically 
minimize the contributions from bad periods both for the 
coincidence search and for the upper limit evaluation. 

Before we start describing the procedure two comments 
are in order : 

1) when segmenting the data in subperiods to be treated 
separately, or to run separate optimizations of the search 
parameters, we must ensure that each subperiod be long 
enough to provide a sufficient statistics, and in particular 
to avoid that any particular outlier or temporary noise 
affect in a sizable way the final choices. We found that a 
few days of data is somehow the minimum duration for 
this purpose. 

2) if the procedure is properly devised, the addition of 
any information or data set, however poor its quality 
with respect to the rest of the data, should not reduce the 
quality of the total result. On the contrary, a correct way 
of putting together all the information can only produce 
a better result. 



Subperiod 






Days of 


V< hi > ■ 10" iy 




Ti 


Begin 


- End 


good data 


Explorer 


Nautilus 


#1 


2007 A 


Apr. 16 


- Dec. 5 


162 


4.57 


3.47 


#2 


2007 B 


Dec. 6 - 


Dec. 31 


12 


4.21 


5.04 


#3 


2008 


Jan. 1 - 


Dec. 31 


232 


4.43 


4.36 


#4 


2009 


Jan. 1 - 


Dec. 31 


242 


3.82 


2.52 


#5 


2010 


Jan. 1 - 


Jun. 10 


113 


3.98 


2.39 



TABLE I: The time stretches Ti{i = 1...5) used to run dif- 
ferent optimization in our analysis. < h'$(t) > is the long 
term average of the amplitude noise 



A. Data selection 

The data were selected with different cuts, applied 
both to the data stream and to the list of outliers. All 



criteria, studied in the past and in use for several years, 
were a priori chosen and blindly applied. The vetos that 
cause elimination of entire periods of data stream include: 

• when acquisition flags or operator's notes are 
recorded, indicating bad or suspected periods (e.g., 
cryogenic refills, activity around the detector....) 

• when the noise of the filtered data, averaged over 
10 minutes, rises above a given value (about 5 times 
the long term average) 

• when the reference tone, a monochromatic signal 
monitoring the gain of the electronic chain, falls 
outside a given range 

• when an excessive amount of wide-band noise is 
present. Wide- band noise, usually of electronics 
origin, is monitored on two frequency bands, above 
and below the useful bandwidth of the detector. 

• when auxiliary channels exhibit mean values above 
predetermined levels. Auxiliary (or veto) channels 
include seismic monitors, SQUID locking working 
point, nitrogen (on Explorer) and helium flow and 
more. 

These cuts reduce the amount of available data for the 
coincidence analysis to 761 days, i.e. two thirds of the 
1152 days of total observation period. The main con- 
tribution to these cuts is due to operations of cryogenic 
maintenance (liquid helium refills) that we chose to per- 
form in different times on the two detectors, so that at 
least one were always operational. 

On these data, an automatic event finder procedure 
selects the "outliers". All data points remaining above a 
chosen threshold are grouped in one event. An event can 
extend over more than one group if the signal falls below 
threshold for a time shorter than the dead time, set to 
1 s. For the class of short signals discussed above, the 
shape of the event is mostly due to the antenna response 
function (see e.g. figH]). Each event is then character- 
ized by the time and amplitude of the largest sample. 
Further characteristic are recorded for each event, such 
as: starting time, total time length, integrated amplitude 
of the samples above threshold, average noise before the 
event. We set for this selection a threshold at critical 
ratio CR=5 with respect to the average noise level, con- 
tinuously updated. 

A further selection was then applied to the outliers, in 
order to implement other cuts: 

• an event should remain above threshold for a time 
consistent with its amplitude (the decay time for 
the filtered data is the inverse of the detector band- 
width). 

• cosmic ray showers are known to produce short 
bursts of excitation in the antennas. The events 
must not be in coincidence with a shower, as 
recorded by detectors installed above and below 
both antennas [III [l7| ■ 
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Subperiod 


Duration 


AT 


AT 




(day) 


in 10 4 shifts 


/ [day ■ shift) 


#1 : 2007 A 


162 


5,635,671 


3.48 


#2 : 2007 B 


12 


640,143 


5.33 


#3 : 2008 


232 


6,667,062 


2.87 


#4 : 2009 


242 


3,061,314 


1.26 


#5 : 2010 


113 


1,042,858 


0.92 



TABLE II: Accidental (time-shifted) coincidences (N acc ) in 
the 5 subperiods analyzed, obtained with the lowest threshold, 
hrss = 3.56 ■ 10~ 20 s 1//2 , used in this analysis. 



These two selections veto a very small fraction of the 
events, usually less than 0.1%. 



III. DETECTORS CHARACTERIZATION 

In order to perform a "fair" search for coincidence, all 
"human handles", i.e. adjustable parameters, must be 
a priori set before starting the search. To this purpose, 
we have applied a very large number of software injected 
events to determine the efficiency of both detectors to 
short bursts of GW. Likewise, the background of acci- 
dental coincidences is determined via a large number of 
time shifts. For a given level of accidentals, a priori set, 
the efficiency is then maximized with a proper choice of 
the thresholds. Only at this point, we can "open the 
box" , i.e. look at the zero delay coincidences and as- 
sess its significance. The following subsections give some 
details about this characterization procedure. 



A. Accidentals 

The evaluation of the expected background of acciden- 
tal coincidence was performed with the usual method of 
the time shifts. The lists of events, one for each detector, 
extracted from the data as described in sectHIl were com- 
pared after shifting the time stamp of one of them. The 
event times of one detector were delayed, with respect to 
the other ones, 10,000 times in steps of 1.5 seconds (i.e. 
between ±7, 500s, excluding the zero time shift). This 
value is larger than the dead time [2(| inserted by the 
event finder. The search of coincidences in each of the 
10,000 cases, performed with a time window of 15 ms as 
discussed in sect llll Cl produces the data base of unphys- 
ical coincidences from which we learned the background 
characteristics. Tab. II summarizes the results. 



B. Software injections 

Large sets of software injections were performed in or- 
der to determine the efficiency of the detectors to delta- 
like signals of different amplitudes. As mentioned above, 
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FIG. 2: Excitation of the detectors, as seen in the filtered 
data, due to two large cosmic-ray showers: an event of Ex- 
plorer 2008 and an event of Nautilus 2007. Only the max- 
imum value of these signals can be interpreted in terms of 
the assumed /i r3S excitation; the shape of the pulse is due 
to the antenna response. The two plots are so dissimilar be- 
cause of the differences in the bandwidths (see fig[l]): while 
Nautilus has its sensitivity around one main frequency, Ex- 
plorer is most sensitive on two frequencies and therefore its 
time-domain response exhibits beats. 



the extensive cosmic ray showers excite the bars, closely 
approaching the effect of a short GW burst. We took 
advantage of this feature and used real signals, observed 
in coincidence with some particularly intense cosmic ray 
shower, as the prototype signal to be used for software 
injections (see figH] for an example of the signals ap- 
plied). These signals, actually oversampled at 50 kHz, 
were scaled to the appropriate values of amplitude and 
added to the filtered data of each detector. This tech- 
nique is much faster than that generally used, where 
one first generates the h-reconstructed data stream, then 
adds the injections to this stream and finally re-filters 
and searches for the events. We validated our method by 
applying both techniques to a one-day sample of data, 
finding a very good agreement. The times of the injec- 
tions were pseudo random, because we avoided injections 
too close to the beginning or end of each period of good 
data, and required a minimum distance of 10 seconds 
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FIG. 3: Efficiencies of Explorer and Nautilus in 
2009. The four lines refer to injections with h r33 — 
(7.97, 9.76, 11.27, 12.6) ■ 10~ 20 s 1/2 . 

between two adjacent injections. Moreover, we added a 
delay, randomly chosen in the [—2.3, 2.3]ms interval to 
the injection time of Explorer, to simulate the time of 
flight of a possible GW signal of unknown direction. 

We injected signals of 10 different amplitudes, in the 
range of h rss [7.97 ~ 12.6] ■ 10~ 20 s 1/2 , at a rate of about 
90 injections per day. 

C. Efficiency 

The usual event finder routine was then applied to the 
data containing the injected signals. For each sub-period 
and for each level of injected signal, efficiency charts were 
produced, displaying the percentage of detected signals 
with amplitude exceeding any given value. Figj3] is a 
sample of such charts, showing the efficiencies in the sub- 
period 2009 for each antenna. 

The injections also allow us to determine the time re- 
sponse of the detectors, and guided us in choosing the 
best coincidence time window to be applied. We found 
that a coincidence window of ±15 ms assures an efficiency 
very close to 1 for delta-like signals, even at the lowest 
injected amplitude: indeed, the measured probability of 



missing a coincident event with the chosen window of 15 
ms is less than 1 • 10~ 4 . Besides, the chosen window is 
sufficiently wide to also accommodate, without signifi- 
cant losses of efficiency, other classes of signals [2fJ for 
which the detectors time response might be different. 



D. Receiver Operating Characteristics 

Efficiency and accidentals vs threshold amplitude com- 
pletely characterize a detector. These two classes of in- 
formation can be summarized in the Receiver Operating 
Characteristics or ROC. 

It is worthwhile, in view of what follows, to briefly re- 
call the procedure to generate a ROC: for each injected 
amplitude we sweep the threshold amplitude and we look 
up both the efficiency and the rate (or the total number) 
of accidentals. By eliminating the threshold value be- 
tween these two relations, we derive a curve [efficiency 
vs event rate], that constitutes the ROC for that given 
signal amplitude. 

In figUwe show an example of ROCs, for both Explorer 
and Nautilus, relative to year 2010. 

It is to be remarked that, despite the fact that the 
detector hardware was virtually unmodified in all sub- 
periods, the ROCs do vary, especially for Nautilus, from 
one subperiod to another (see figE]). 

As we are interested in the operation of both antennas 
as one detector, we can extend the concept of ROC to a 
coincidence search. In this case, in order to vary efficiency 
or accidental rate, we can act on either threshold, so that 
there exists an infinity of threshold pairs that can provide 
the same characteristics: we could have therefore infinite 
ROC curves for the same signal amplitude. However, 
keeping in mind that our aim is to maximize the efficiency 
of detection for a given accidental rate, we can focus our 
search on finding the pair of thresholds that gives the 
best efficiency for each value of accidentals. The ROCs 
for the Explorer + Nautilus compound observatory are 
therefore obtained with the following procedure: 

• we choose a set of M threshold values and we sift 
through our data with a matrix of M*M thresholds: 
both in the list of events found with the injections 
where we demand triple coincidences (tsx , tjsiA and 
ti n j), and in the set of shifted coincidences, for the 
accidentals. 

• in this way we create a M*M matrix with values of 
efficiency and accidentals for each threshold pair. 

• for each of N chosen values of accidentals, we search 
the matrix for those intervals that contain that 
value of accidentals. We interpolate in those inter- 
vals to find the value of thresholds and efficiency. 

• finally, we compare these values and choose that 
with the largest efficiency 
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FIG. 6: ROCs for the combined observatory (Explorer + Nautilus). Left: in 2010, at the four injection amplitudes considered 
in figlU right: at an injected amplitude h raa = 9.76 ■ 10 _20 s 1//2 in the five subperiods. 



In our search, we used N=M= 100 and repeated the pro- 
cedure for the 5 subperiods and for each of the 10 values 
of injected amplitudes. 

Some of the ROCs for the observatory, obtained with 
this procedure, are shown in figlU Each point of these 
curves represents the threshold pair that produces the 
desired value of accidentals with the best possible effi- 
ciency. This procedure yields approximated values for 
the data points, as they are obtained via interpolation: 
for this reason, the search was later refined around the 
selected threshold values. 



IV. COINCIDENCE SEARCH 

In order to perform the true-time (on-source) search, 
we must decide upon a unique set of 5*2 thresholds to 
be applied to the 2 detectors in each of the 5 subperi- 
ods. This set must provide the desired number of total 
accidentals (0.1) while achieving the maximum possible 
efficiency for GW signals. 

The ROCs for the coincidences, previously determined, 
specify, in each subperiod and for each injection am- 
plitude, what are the thresholds capable of obtaining a 
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FIG. 7: Maximum efficiency e achievable with a background 
set at 0.1 accidentals in the entire observation time. The top 
line shows the efficiency with thresholds optimized at each 
abscissa point (injected amplitude), while the other curves are 
produced with different choices of a unique set of thresholds. 

given value of accidentals with the maximum possible ef- 
ficiency. Next step is to find how to distribute the total 
number of accidentals between the different periods in 
order to maximize the total efficiency defined as 

e = £^P (1) 

j=i 

where £j, are the efficiency and duration of the subpc- 
riods and T = £\ Tj is the total observation time. 

We remark that this procedure pins down a different 
set of thresholds for each considered injection amplitude. 
In figlZl the data points show the results of this opti- 
mization: each data point is obtained with its own op- 
timized set of 10 thresholds; we call this curve "com- 
posite efficiencies". However, as the coincidence search 
has to be performed only once, we need a strategy to 
select a unique set of thresholds. FigJT] also shows three 
curves describing the efficiency at all amplitudes for three 
selected sets of thresholds, namely those optimized for 
(7.97, 9.43, 11.96) -10~ 20 s 1/2 . We make here no assump- 
tion on the amplitude distribution of the GW signals we 
search for; therefore we selected a threshold set that best 
approaches the curve of "composite efficiencies" at all 
amplitudes, and in particular at the smaller ones (that 
are, in such a search, the most probable). 

Clearly, the set of thresholds chosen for h rss = 9.43 ■ 
10 _20 s 1 / 2 (values detailed in table III) is the one that 
performs best and was therefore selected. 

Table IIVI shows how the overall background was dis- 
tributed and how the efficiency of detection at several 
signal amplitudes changed over the 5 subperiods of the 
search. We note that the optimization procedure auto- 
matically weights the subperiods according to the data 



Period 




-30,1/2] 






Expl 


Naut 


#1 


2007 A 


9.00 


8.52 


#2 


2007 B 


12.3 


9.60 


#3 


2008 


10.8 


19.3 


#4 


2009 


8.12 


8.17 


#5 


2010 


8.19 


8.29 



TABLE III: The set of thresholds chosen for the final coinci- 
dence search: these values are optimal for an injected signal 
with h raa = 9.43 ■ 10 _20 s 1//2 , but well approach the best pos- 
sible efficiency at all signal amplitudes, as shown in figO 







Su 


b Period 






2007 A 


2007 B 


2008 


2009 


2010 


Accidentals 


0.0206 


0.0 


0.0037 


0.0585 


0.0172 


Injected h rS s 
[s 1 / 2 ] ■ lO" 20 


Efficiency 


7.97 


0.0935 


0.0 


0.0 


0.2066 


0.1520 


8.36 


0.1826 


0.0018 


0.0 


0.3787 


0.3162 


8.73 


0.2936 


0.0018 


0.0 


0.5542 


0.5137 


9.09 


0.4167 


0.0036 


0.0 


0.7063 


0.6856 


9.43 


0.5438 


0.0100 


0.0 


0.8143 


0.8081 


9.76 


0.6528 


0.0190 


0.0 


0.8840 


0.8869 


10.54 


0.8476 


0.0805 


0.0 


0.9533 


0.9712 


11.27 


0.9358 


0.2081 


0.0 


0.9758 


0.9920 


11.95 


0.9709 


0.4009 


0.0 


0.9843 


0.9963 


12.60 


0.9834 


0.6091 


0.0 


0.9890 


0.9980 



TABLE IV: Accidentals and efficiencies at various amplitudes 
with the chosen set of thresholds (see tab. III). Note that 
2007B gives no contribution to the accidental background, 
while for 2008 the efficiency results set to zero at all ampli- 
tudes. 



quality, virtually "turning off", without any manual ad- 
justment, the noisiest periods, i.e. 2007B and 2008: 
where we have a noisier detector, there we get little or 
no contribution to the coincidence search. 

When we finally applied this set of thresholds to the 
on-time data, no coincident events were found, thus re- 
turning a null result. 



V. UPPER LIMITS 
A. Method 

We now describe the procedure employed to compute 
the upper limits on the rate of incoming GW short bursts 
for a set of possible signal amplitudes: this procedure is 
separately applied to each of the 5 subperiods in which 
the entire observation time T was segmented. These re- 
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suits are then combined and an overall 95 % bayesian 
upper limit is determined at each signal amplitude. 

We remark that, when we compute the upper limit 
(UL) for a given GW amplitude, we are assuming the hy- 
pothesis that only signals of that very amplitude could 
reach the Earth. This means that each point in a UL 
curve is independent of any other point, and its determi- 
nation can be independently optimized. 

The handles we have for this optimization are, just 
as in the coincidence search previously described, the 
thresholds to be applied to the data: varying the thresh- 
olds allows us, in turn, to change: 

• the background, i.e. the rate ro, or the mean 
total number /io = tq ■ T, of accidental coinci- 
dences. These are, as before, estimated with the 
time-shifted data fsect llll X)) . 

• the efficiencies e, as computed with the software 
injections fsect llll Cj) . 

The output of this search is the estimated maximum rate 
r of incoming GW signals, at any signal amplitude h rss 
or, equivalently, the total number of detected GW signals 
fi = e ■ r ■ T. The optimization consists in choosing the 
thresholds potentially capable of producing the best, i.e. 
lowest, upper limit r{h rss ). 

We note that the optimization procedure is different 
from that employed in sect II VI for the coincidence search: 
in that case we looked, at each amplitude, for the thresh- 
olds that would yield the best efficiency for a given (0.1 
events) background of accidentals. Here, not being tied 
to a pre-fixed value of accidentals, we can choose the pair 
efficiency-background that optimizes our result. 



B. The relative belief updating ratio R 

The quantity we need to compute and optimize, for 
each value of assumed signal amplitude, is the relative 
belief updating ratio R, i.e. the ratio of the likelihood 
P (Mo + I 1 ! N) that the N coincident events found be due 
to the presence of a given number /i of GW events, to the 
likelihood P(fio,N) of a mere accidental background. 

By assuming, as usual, that the number N of coinci- 
dences found obeys the Poisson statistics, we then can 
write the likelihood in the presence of a rate r (corre- 
sponding to a detectable number fj) of GW events as: 

P( M o+M,iV)= (Mo+/ y (2) 

The same relation, with \x = 0, describes the likeli- 
hood P(fj,o,N) of mere background. Therefore, the rela- 
tive belief updating ratio can be written, in terms of our 
parameters, as: 



R(r) 



Mo" 



sr 
'•() 



N 



The determination of the maximum rate r(h rss ) of 
GW, reaching the Earth with a given amplitude, requires 
the elaborate procedure outlined below. For sake of clar- 
ity, we describe it for a fixed value of h rss , implying that 
it is then repeated for all amplitudes of interest. 

Requiring a 95% confidence limit means finding the 
particular value r* such that R(r*) = 0.05. Eq. [3] shows 
that R also depends on other parameters, so that an 
optimization is possible: we can vary our handles (i.e. 
the thresholds) till we find the minimum value among all 
r* . The functional dependence of R on the thresholds is 
due to two competing effects: by raising the thresholds, 
we decrease both r (that decreases R) and e, i.e. \i (that 
increases R). The result can't be analytically predicted, 
and a numerical search should be done by varying the 
thresholds. Actually, part of this work has already been 
done in computing the ROCs (see sect HIIDj) : we found 
there the optimal efficiency for each number (or rate) 
of accidentals. This simplifies our search: rather than 
probing the entire (ro, e) plane, we only need to compute 
R(r|ro, s) along the ROC curve corresponding to the h rss 
considered, as in fig [SI We now have a family of curves 
i?(r) depending on one parameter: the (ro, e) pair. Of all 
these curves, we select the one where the relation R(r*) = 
0.05 is achieved with the lowest value r* . 

One last ingredient is missing for this calculation: the 
total number TV of coincidences: it is provided by the " on- 
time" analysis that cannot be performed before setting 
all the search parameters. We have then implemented 
the following workaround: for each pair (/zo,e) from the 
ROC. Wc solve for r* the relation R(r*) = 0.05 at all 
possible values of iV acc . We then compute the weighted 
average of these r's: 



(4) 



In practice, the sum is truncated when 

M 

J2P(Vo,N acc )>l~l-10- 10 (5) 



N acc 



e~ re e~ T 



(3) 



Finally, we choose the threshold pair that yields the 
minimum f*, i.e. that minimizes the expected UL based 
on our efficiency and accidental rate. The search is then 
repeated at a different value of GW amplitude, until the 
curve r(h rss ) is traced. 

The entire calculation is an " a priori" procedure, per- 
formed without any knowledge of the on time coinci- 
dences. 



C. Upper Limit evaluation 

The above procedure was applied to find the optimal 
threshold pairs in the 5 time subperiods and for 19 dif- 
ferent injection amplitudes [25| Only at this point, we 
had the right to "open the box" and find the on time 
coincidences. 
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FIG. 8: Thresholds of Explorer and Nautilus resulting from the UL optimization. 
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FIG. 9: .R-curves at h ras = 5.63 ■ IQ~™ (a) and 7.97 ■ lO - * (b) s 1/2 in the five sub-periods of data taking 
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FIG. 10: Total .R-curves at different values of injected h r 



Figure 8 shows the thresholds for Explorer and Nau- 
tilus resulting from the UL optimization. 

Figures 9a and 9b show, as an example, the -R(r) curves 
at h™i = 5.63- 10" 20 and 7.97- lO^s 1 / 2 in the five subpe- 
riods, as well as the total R(r), obtained by multiplying 
the curves of all subperiods. In fig. 10a and 10b we report 
the total i2(r) computed at various h™J a . 

The standard procedure, at this point, calls for eval- 
uation of the UL as the product of R(j) and a prior, 
containing all our previous knowledge on ULs. The best 



prior is, in principle, a combination of all previously com- 
puted ULs (e.g. 1, i). However, this cannot be used, 
due to different meanings and methods of these ULs, as 
discussed in the next section. We are then left with the 
choice of a purely theoretical prior, as we make no as- 
sumptions on the source location, polarization, sky dis- 
tribution etc., it is reasonable to assume a flat prior. Our 
95% UL is then simply obtained by picking the value of 
r that yields #=0.05. 

Figure [TT] shows the on time results for the UL at dif- 
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FIG. 11: Comparison between the upper limit computed with 
the optimized procedure described in the text (opt) and with 
a procedure with fixed thresholds, those defined by the co- 
incidence search (fix). There is an evident large gain at low 
amplitude. 
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FIG. 12: Comparison between our 95% Upper Limit and pre- 
viously published results. All LIGO and Virgo results, and in 
particular the solid curve (LV2), refer to 90% ULs. 



ferent values of h rss ■ The statistical uncertainties on the 
determination of efficiencies and accidentals (fig) causes 
on the evaluation of the ULs a relative error at 1 a rang- 
ing from 4% to 0.4% for h rss going from the lowest to the 
highest value, respectively. The various parameters en- 
tering the final evaluation of the UL, namely the number 
of on-time coincidences iV co j, the estimated accidentals 
Ho and the computed efficiences, are reported in tab.V 
for all the values of injected h rss and all five subperiods. 



VI. DISCUSSION AND CONCLUSIONS 



A. Comparison with other experiments 



FigHU compares our present results with some other 
published in the past: the curve labeled "IGEC 1" is 
from @, "IGEC 2" is from [|, "S2" is the UL for 1 ms 
gaussian pulses from the LIGO S2 run|2l| , " S4" ," S5" and 
"LV2" are for Q=9 sin-gaussian pulses at 1053 Hz from 
|3, H and 0. 



In the present calculation of UL there are several 
choices that make it difficult to compare with previous 
results. 

If we consider the above cited upper limits released by 
the Ligo Scientific Collaboration (LSC), we can see that 
no specific optimization was carried out for the UL: all 
the analysis parameters, and in particular the thresholds, 
were set for the coincidence search. The number of on- 
time coincidences found was directly used to compute 
the UL for that run. We show in fig [TT] what the result of 
such an analysis strategy would be on our data: the curve 
labeled "opt" represents the optimized UL computed in 
the previous section while the curve labeled "fix" is the 
result we would obtain with the thresholds determined in 
the coincidence search. The improvement in sensitivity, 
especially at low amplitudes can be clearly seen. 

A more relevant issue regards the meaning of the vari- 
able h rss , i-e. the abscissa of the UL plot. In the two 
IGEC searches 0, H[ the efficiency of the detectors was 
not considered, and the UL was plotted vs the amplitude 
threshold used in the coincidence search. What the curve 
really meant was then the UL on GW rates detectable 
by the observatory with that threshold, rather than the 
incoming rate. For a fair comparison, at least the effi- 
ciency of the detectors should be considered: for a given 
amplitude ho, when the threshold is set at that very am- 
plitude ho, the efficiency is roughly (l/2) n where n is the 
number of detectors. The IGEC UL values should there- 
fore be increased by at least a factor 4 to convert to an 
incoming rate. 

On the other hand, the LSC includes in its analysis a 
model describing the distributions of the GW incoming 
signals, by folding into the calculation a factor f(0, 4>) < 1 
to account for isotropic direction and random polariza- 
tion of the wave. The h rss in LIGO's UL is then the 
maximum amplitude detectable by the observatory. Be- 
sides, they computed a 90% UL, while ours is at 95%: 
according to eq.4.1 of 3, in order to convert their ULs 
at 95%, they should be increased by a factor ~ 1.3. It's 
clear that the combination of these two differences, one 
of which would lead to a decrease of the UL, the other 
one to an increase, would anyhow not change much the 
comparison with our results. 
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2007 A 




2007 B 
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hrss 


N coi 


Ho 


e 


N coi 


Ho 


e 


N coi 


Ho 


e 


N coi 


Ho 


e 


N coi 


Ho 


e 


2.52 


565 


553.9 


.0245 


66 


63.04 


.0179 


639 


663.2 


.0159 


283 


303.4 


.0193 


94 


103.2 


.0161 


3.56 


565 


553.9 


.1396 


66 


63.04 


.0868 


638 


663.2 


.0910 


273 


295.5 


.1909 


94 


103.2 


.2227 


4.36 


95 


85.67 


.1432 


66 


63.04 


.2139 


638 


663.2 


.2272 


109 


130.8 


.3536 


40 


45.68 


.4378 


5.04 


57 


48.94 


.2532 


32 


23.92 


.2612 


638 


663.2 


.3784 


30 


37.01 


.3699 


18 


24.30 


.6078 


5.63 


20 


18.74 


.3090 


26 


21.31 


.3864 


390 


408.2 


.4192 


15 


15.95 


.4611 


5 


3.856 


.5261 


6.17 


5 


7.177 


.3600 


25 


19.89 


.5002 


111 


97.72 


.3133 


1 


3.208 


.4180 


2 


1.536 


.6222 


6.67 


2 


3.588 


.4339 


7 


6.881 


.4406 


46 


41.23 


.3354 


1 


1.457 


.5283 


1 


.7959 


.7130 


7.13 


1 


2.162 


.5288 


5 


4.439 


.5068 


30 


25.37 


.4041 





.9320 


.6522 


1 


.5365 


.8032 


7.56 





.7643 


.5472 


3 


2.930 


.5518 


19 


12.93 


.4341 





.5226 


.7189 





.3088 


.8393 


7.97 





.5404 


.6346 


2 


1.890 


.5836 


8 


8.396 


.4946 





.3342 


.7762 





.2028 


.8728 


8.36 





.3344 


.6941 


1 


1.935 


.6806 


7 


6.411 


.5696 





.2437 


.8347 





.1403 


.9019 


8.73 





.2242 


.7465 


1 


1.646 


.7292 


3 


4.395 


.6091 





.1643 


.8606 





.0749 


.9031 


9.09 





.1503 


.7892 


1 


1.401 


.7623 


2 


3.095 


.6403 





.1297 


.8899 





.0589 


.9289 


9.43 





.1063 


.8281 


1 


1.369 


.8166 





2.122 


.6620 





.1079 


.9126 





.0477 


.9497 


9.76 





.0712 


.8526 


1 


1.248 


.8491 





1.805 


.7164 





.0787 


.9217 





.0387 


.9650 


10.54 





.0491 


.9257 





.2660 


.7679 





1.024 


.7810 





.0545 


.9548 





.0223 


.9834 


11.27 





.0296 


.9537 





.1186 


.8439 





.5361 


.7897 





.0387 


.9652 





.0135 


.9892 


11.95 





.0140 


.9647 





.0594 


.8952 





.5361 


.8967 





.0305 


.9755 





.0084 


.9919 


12.60 





.0113 


.9790 





.0265 


.9210 





.5361 


.9455 





.0254 


.9821 





.0068 


.9959 



TABLE V: Characteristic parameters computed in order to evaluate the R(r) curves in the 5 subperiods: N CO i is the on time 
number of coincidences detected, is the estimated average background and e is the efficiency of detection at each particular 
value of h rss 



B. Does the sub-periods segmentation pay off ? 



C. Conclusions 



It is reasonable to question whether a unique search 
over the entire observation period T would yield a sim- 
ilar or better result with respect to our choice of seg- 
menting the analysis in 5 sub-periods. In other terms, 
how does the global Rg(r) compare vs the product 
R-t(0 = n^>( r ) °f 5 separate update ratios? An ex- 
act answer can be given in the simple case where we 
segment T in two subperiods having the same character- 
istics, namely accidentals ho and efficiency e, assumed 
constant. In this case, the global (one period) R is given 
by eqUJ while the product of the two R's for the subpe- 
riods Ti,T 2 is: 



R 



w = n 



{Ht + Ho,t) 1 e 



(0) 



As Ho — H0.1+H0.2 (and same for /i) and N — N1+N2, 
one can expand eq[6l and prove Rt = Rg , for any choice 
of Ti,T 2 , iVi,JV 2 . 

For the more general case of two non homogeneous 
subperiods, although we lack an algebraic proof, exten- 
sive numerical investigation has shown that we should 
always expect Rt <Rg- 



In this paper we analyzed 3 years of almost continuous 
data from the two resonant gravitational wave detectors 
Explorer and Nautilus. The period examined spans from 
the end of the IGEC2 four-detector analysis to the de- 
commissioning of Explorer. Both the search for coinci- 
dences with low false alarm rate and the evaluation of 
the upper limit have been performed employing a novel 
type of analysis, with optimization of the thresholds of 
each detector separately for each intermediate task. This 
method has proven successful in obtaining better results 
(see for instance figfTT]) as well as for handling non sta- 
tionarities in the detectors behavior. As an example, we 
recall the noisy period of 2007B: in a search with the 
usual procedure, that period would be discarded, or its 
large number of events would negatively affect the statis- 
tics of the remaining, better data. In our case, the op- 
timized procedure automatically takes care of the higher 
noise and reduces the weight of that period on the final 
results. Indeed, its contribution to both the coincidence 
search and the upper limit evaluation is hardly notice- 
able. 

The upper limit computed on the basis of our data 
cannot compete with those of the more sensitive inter- 
ferometric detectors, that extend down to much smaller 
Nevertheless, the length of our data col- 
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lection lead us to expect that we could improve upon 
the UL set by LIGO S5, at amplitudes of the order of 
hras ~ lO-^s 1 / 2 : indeed we did obtain a better UL than 
other previously published. However, while in the process 
of analyzing our data, a new, improved UL was released 
by the LSC- VIRGO collaboration: combining the data 
of the S5/VSR1 and S6/VSR2-3 runs, the extended data 
taking allowed them to set a better limit also at higher 
amplitudes. 

The procedure detailed here could be profitably used 
in future searches, where better sensitivity of the detec- 
tors would yield even more significant ULs. Infact, we 
demonstrated (see fig fTTj) that this procedure grants a 
substantial improvement in the evaluation of the Upper 
Limit, up to two orders of magnitude at low amplitudes, 
with respect to the standard way of computing it. 



Acknowledgments 

Explorer ceased taking data on June 11th, 2010, con- 
cluding an activity that spanned over a quarter of a cen- 
tury. We take the opportunity of this last scientific paper 
produced with its data to thank CERN as an institution 
for hospitality and provision of facilities, first and fore- 
most the cryogenic fluids. We also thank the numerous 
CERN officials and staff that provided, over such a long 
period, with logistics and with daily as well as emergency 
help. 

We are grateful to the technicians (too numerous to name 
them all) that, over the years, helped us maintaining both 
detectors in operations and, to these days, still work to 
keep Nautilus on the air. 



References 

[1] Z.A.Allen et al., Phys.Rev.Lett. 85, 5046 (2000). 
[2] P.Astone et al., Phys.Rev.D 68, 022001 (2003). 
[3] P.Astone et al., Phys. Rev. D 76, 102001 (2007). 
[4] See for instance J.Abadie et al., Phys.Rev.D 85, 082002 
(2012). 

[5] E.Amaldi et al., Nuovo Cim.C 7, 338-354 (1984). 

[6] P.Astone et al., Phys.Rev.D 47, 362-375 (1993). 

[7] P.Astone et al., Phys.Rev.Lett. 91, 111101 (2003). 

[8] P.Astone et al., Europhys.Lett. 16, 231-235 (1991). 

[9] P.Astone et al., Astropart. Phys. 7, 231-243 (1997). 
[10] P.Astone et al., Class.Quant.Grav. 19, 5449-5463 (2002). 
[11] P.Astone et al., Class.Quant.Grav. 23, S169 (2006). 
[12] P.Astone et al., Phys.Lett.B 385, 421-424 (1996). 
[13] J.Zendri et al., Class.Quant.Grav. 19, S1925 (2002). 
[14] M.P.McHugh et al., Class.Quant.Grav. 22, S965 (2005). 
[15] E.Coccia et al., Nucl.Inst.Meth.Phys.Res.A 355, 624-631 
(1995). 

[16] P.Astone et al., Phys.Rev.Lett. 84, 14 (2000). 



[17] P.Astone et al., Astropart. Phys. 30, 200 (2008) . 
[18] P.Astone et al., Phys.Lett.B 499, 16-22 (2001). 
[19] M.Bassan et al., Nucl.Inst.Meth.Phys.Res.A 659, 289- 
298 (2011). 

[20] P.Astone et al., Phys.Rev.D 82, 022003 (2010). 

[21] B.Abbott et al., Phys.Rev.D 72, 062001 (2005). 

[22] B.Abbott et al., Class.Quant.Grav. 24, 5343 (2007). 

[23] B.Abbott et al., Phys.Rev.D 80, 102002 (2009). 

[24] Traditionally, the GW resonant detectors community 
indicates the antenna sensitivity, i.e. the long term 
quadratic average of the filtered output hs (t) , in terms of 
a "detection noise temperature" T e ff. These two quan- 
tities are related, for our detectors, by < hg >= 6.35 • 
10 — 35 T e ff. In these units, the average noise levels of 
table I spanned between 0.9 and 4 mK. 

[25] Beside the 10 values of injected amplitude used in the co- 
incidence search, we also explored the small signal region 
repeating the injection at 9 additional amplitudes in the 
range [2.52 ~ 7.56] ■ 10" 2 V /2 . 



